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ABSTRACT 

We present new sub-arcsecond (0.7") Combined Array for Research in Millimeter-wave As- 
tronomy (CARMA) observations of the 1.3 mm continuum emission from circumstellar disks 
around 11 low and intermediate mass pre-main sequence stars. High resolution observations for 
3 additional sources were obtained from literature. In all cases the disk emission is spatially 
resolved. We adopt a self consistent accretion disk model based on the similarity solution for 
the disk surface density and constrain the dust radial density distribution on spatial scales of 
about 40 AU. Disk surface densities appear to be correlated with the stellar ages where the 
characteristic disk radius increases from 20 AU to 100 AU over about 5 Myr. This disk expan- 
sion is accompanied by a decrease in the mass accretion rate, suggesting that our sample disks 
form an evolutionary sequence. Interpreting our results in terms of the temporal evolution of 
a viscous a-disk, we estimate (i) that at the beginning of the disk evolution about 60% of the 
circumstellar material was located inside radii of 25-40 AU, (ii) that disks formed with masses 
from 0.05 to 0.4 M© and (iii) that the viscous timescale at the disk initial radius is about 0.1-0.3 
Myr. Viscous disk models tightly link the surface density with the radial profile of the disk 

viscosity v(R) cx B? . We find values of 7 ranging from -0.8 to 0.8, suggesting that the viscosity 
dependence on the orbital radius can be very different in the observed disks. Adopting the a 
parameterization for the viscosity, we argue that a must decrease with the orbital radius and 
that it may vary between 0.5 and 10~ 4 . From the inferred disk initial radii we derive specific 
angular momenta, j, for parent cores of (0.8 — 4) x 10 -4 km/s pc. Comparison with the values 
of j in dense cores suggests that about 10% of core angular momentum and 30% of the core 
mass are conserved in the formation of the star/disk system. We demonstrate that the similarity 
solution for the surface density for 7 < can explain the properties of some "transitional disks" 
without requiring discontinuities in the disk surface density. In the case of LkCa 15, a smooth 
distribution of material from few stellar radii to about 240 AU can produce both the observed 
SED and the spatially resolved continuum emission at millimeter wavelengths. Finally we show 
that among the observed sample, TW Hya is the only object that has a disk radius comparable 
with the early solar nebula. 



Subject headings: 

1. Introduction 

Spatially unresolved observations of the in- 
frared and mm-wave emission from nearby pre- 
main sequence stars surrounded by disks suggest 
that most of the circumstellar dust dissipates 
on timescales between 1 and 10 Myr (see, e.g., 



Hernandez et al. 1120071 ). Nevertheless, it remains 



very uncertain how disk evolution proceeds in in- 



dividual systems and, in particular, whether all 
circumstellar disks give rise to planetary systems. 
Over the last ten years a large number of circum- 
stellar disks in nearby star forming regions have 
been observed using long baseline millimeter and 
sub-millimeter interferometers. These observa- 
tions have spatially resolved the disk emission to 
infer the radial distribution of gas and dust. How- 
ever, at the distance of the nearby star forming re- 
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gions (100- 200 pc), the typical angular resolution 
of 1.5 "-3" fRodmann et al. Il2006t lAndrews et aT 



20071 : iKitamura et al. 1 120021 : iDutrev et al. I ll996) 



could not constrain the detailed structure of disks, 
which typically have radii of only a few hundred 
AU. 

Higher angular resolution mm-wave observa- 
tions remain challenging and the number of disks 
observed at angular resolution higher than 1" is 
still very small and essentially restricted to more 
massive and luminous pre-main seq uence circum 
stella r disks such as LkHa 
20081 ). HD 163296 (llsella et al 



330 JBrown et al. 



2007|). AB AuT 



MWC 480. DM Tau and Lk Ca 15 (IPietu et al I 



20051 12009 . l2007h . CQ Tau dTesti et al. I l2003h. 
DL Tau, UZ Tau BP Tau and G M Aur (jSimon et ah 
2000h . TW Hva (IWilner et al. l [2000h and GG Tau 



([Guilloteau et al. II1999T ). Even for this small sam- 
ple the radial distribution and kinematics of the 
circumstellar material vary considerably from ob- 
ject to object. Since the observed objects are 
characterized by stellar ages between ~0.1 and 
~10 Myr, which is probably a considerable frac- 
tion of the disk life time, variations in the dust 
properties may also be representative of different 
evolutionary stages. Differences in disk structure 
are believed to result from variation in the total 
angular momenta, masses, chemical compositions 
and magnetic fields, during th e collapse of the 



paren t molecular core (see, e.g.. lHueso fc Guillot 



2005). Disk structure in multiple systems or in 
dense star forming regions can also be influenced 
by the dynamical perturbation induced by close- 
by companions or by strong interstella r radiation 



field (see, e.g, Alexander et al. 20061 ). Detailed 



investigations of disk structure and of the origins 
of any observed difference are clearly necessary 
to improve our understanding of the formation of 
planetary systems. 

Here we present sub-arcsecond observations of 
circumstellar disks around 14 nearby pre-main 
sequence stars. New 1.3 mm continuum obser- 
vations of 11 objects in the Taurus and Ophi- 
uchus star forming regions, CY Tau, DG Tau, 
DM Tau, DN Tau, DR Tau, GO Tau, LkCal5, 
RY Tau, UZ Tau E, GSS 39, SR 24 S, were ob- 
tained with the Combined Array for Research in 
Millimeter Astronom£| (CARMA). For MWC 275 



(HD163296), GM Aur, and TW Hya, we have 
reanalyzed published data from the SMA and 



Plateau de Bure interf erometers ([Hughes et al. 
120081; llsella et al. Il2007l) . For each object we de- 



rive the radial dust distribution by comparing the 
observed dust continuum emission with a self- 
consistent disk model based on the similarity solu- 
tion for the surface density of a vi scous keplerian 



1971 . 



disk (|Lynden-Bell fc Pringle 

In § [5] we summarize the properties of the stel- 
lar sample. The interferometric observations and 
data reduction procedures are described in § [3] In 
§ 2] we discuss the observations. § [5] contains the 
description of the adopted disk model, while the 
results are presented in § [5J The discussion and 
the conclusions follow in § [7] and [8] 

2. The sample 

The 14 stars selected for study are listed in Ta- 
ble [1] together with the adopted luminosities, spec- 
tral types, temperatures and accretion luminosi- 
ties. All are nearby pre-main sequence T Tauri 
stars with known 1.3 mm flux densities in ex- 
cess of 50 mjy, to ensure high signal-to-noise ra- 
tios for the extended disk emission. Two targets, 
UZ Tau E and SR24 S, are members of multi- 
ple systems. The first nine objects of Table [1] 
are located in the Taurus- Auriga star forming re- 
gion, while GSS 39 and SR 24 S are in Ophi- 
uchus. TW Hya is in the homonymous association, 
and MWC 275 (HD 163296) is an isolated Herbig 
Ae star. Based on observations of the molecular 
gas emission, the strong excess continuum emis- 
sion at IR, mm and radio wavelengths appears to 
originate from large amou nts of gas and dust dis- 



Dutrev et al. 1996: Simon et al. 


2000; iQi et al. 


2004; Isella et al. 2007). 



2.1. Stellar properties 

For the objects in Taurus- Auriga and Ophi- 
uchus we assume stellar distances of 140±10 and 
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130±15 pc respectively (see iRebull et al. 1 120041) . 
For TW Hya and MWC 275 we adopt the Hippar- 
cos distances of 56±5 pc (Wichmann et al. 1998) 
and 122±20 pc (van den Ancker et al. 1998) re- 
spectively. 

Stellar ages and masses are derived from 
the H-R diagram using theoretical tracks from 



mg 

D'Antona fc Mazzitellil (|1997l with the 1998 up- 
dated version available on the web; hereafter 
DM97) adopting published spectral types and 
luminosities (see Table [T]). Assuming errors of 
about 30% in the stellar luminosity and half a 
spectral type in the spectral type classification, 
the resulting masses and ages are uncertain by 
30-50%. Although masses and ages are strongly 
dependent on the adopted stellar evolution model 
(Appendix IA"|). the main results of the paper are 
almost independent of this choice, as discusses in 

§13 

Stellar radii based on the effective tempera- 
ture and the bolometric luminosity were combined 
with the derived masses to provide mass accre- 
tion rates M acc , following the relation M acc = 
L aC cR*l '{GM+). Both quantities are listed in Ta- 
ble [TJ together with stellar masses and ages. 

3. Observations 

Interferometric observations of the 11 disks in 
Table[2]were carried out with CARMA, which con- 
sists of six 10 m and nine 6 m antennas, and is lo- 
cated near Big Pine (CA) at an altitude of about 
2200 m. 

The data were obtained between Oct 2007 and 
Apr 2008 using the C and B array configurations 
to provide baseline lengths between 20 and 270 m 
and between 90 and 900 m respectively, corre- 
sponding to angular resolutions of about 0.7" and 
0.4" at 1.3 mm. The CARMA correlator was con- 
figured with two wide bands of 500 MHz each 
and one narrow band of 8 MHz centered at the 
frequency of 230.538 GHz, corresponding to the 
12 CO (2-1) molecular transition. Only continuum 
emission observed in the wide 1 GHz band will 
be discussed here. Each source was observed for 
4-8 hrs to achieve a noise level between 1 and 5 
mjy Source names, array configurations, synthe- 
sized beam sizes, dates and adopted gain calibra- 
tors are summarized in Table [2j 

The data were reduced using the MIRIAD soft- 



ware package. The resulting 1.3 mm maps, ob- 
tained using natural weighting, are shown in Fig- 
ure [U Bandpass calibration relied on observations 
of 3C273 and absolute flux calibration was set by 
observing Uranus. Since we are interested in re- 
solving the disk extended structures we took par- 
ticular care in correcting for the atmospheric tur- 
bulence which may result in artificially extended 
sources. For each set of observations, we cali- 
brated the data by observing every 15 minutes 
a bright unresolved calibrator located within 20° 
from the target. Fitting a bidimensional Gaussian 
profile to the calibrator image we measured the 
seeing and corrected the data using the SEEING 
option of the UVCAL MIRIAD command. Resid- 
ual effects of the seeing are then eliminated by 
rejecting the UV points on which the atmospheric 
turbulence introduces a flux loss higher than 10% 
on the phase calibrator. To check the quality of 
the phase calibration we observed, if available, a 
fainter second point source located within 5° from 
to the science target. We then verified that, after 
applying the calibration solution derived for the 
bright calibrator, we obtain an unresolved image 
of this second calibrator. 

Details of the 1.3 mm SMA dust continuum ob- 
servations and data calibra tion of GM T a u and 
TW Hya can be found in iHughes et al~l (|2008L 



and references therein) . The angular resolution of 
these observations was about 1" at 1.3 mm. Simi- 
larly, details on t he PdB observation s of MWC 275 
can be found in llsella et al. I ((2007). In this case 



the resolution was 2"x0.4". The 1.3 mm maps of 
these three sources are shown in Figure [2] 

4. Observational results 

The observed continuum emission toward each 
of the disks is spatially resolved, centered on the 
position of the parent star and shows an almost 
centro-symmetric surface brightness profile. In 
Table [3] we report the spatially integrated flux 
(column 2), the angular size obtained by fitting 
the continuum map with a bidimensional Gaussian 
profile (column 3), a first order derivation of the 
disk inclination measured from the aspect ratio of 
the Gaussian fitting (column 4), and the position 
angle measured East from North (column 5). Col- 
umn 6 shows the disk radius, R g , defined as the 
radius containing 95% of the observed emission. 
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In practice, R g ~ 1.4 x FWHM measured along 
the apparent disk major axis assuming the stellar 
distances discussed in § O Given the sensitivity of 
the observations and the fact that the dust emis- 
sion is optically thin, this radius provides only a 
rough estimate of the real extent of the disk. 

The measured integrated fluxes are in good 
agreement (20%) with earlier interferometric ob- 
servations at lower angular resolution, suggesting 
that only a small fraction of the flux is emitted 
on completely resolved out disk scales. We there- 
fore assume that the CARMA observations trace 
the bulk of the disk emission. Disk radii range 
from 90 AU (DR Tau) to 320 AU (SR 24 S) and 
the surface brightness slopes vary from steep (e.g., 
DG Tau, DR Tau, DN Tau) to quite shallow (e.g., 
RY Tau, CY Tau, LkCa 15, GSS 39). This latter 
point is illustrated in Figure [21 which shows the 
normalized visibility amplitude V as a function of 
the baseline length between and 270 m for four 
disks. The most extended source, GSS 39 (filled 
circles) , is completely resolved out on a baseline of 
150 m, corresponding to an angular scale of about 
1.8", while the most compact, DR Tau (empty 
squares), is resolved out only on a baseline of 
about 400 m, corresponding to an angular scale 
of -0.7". RY Tau (empty circles) and LkCa 15 
(filled squares) are intermediate cases, resolved on 
angular scales of 0.8" and and 1" respectively. 

Given the high S/N ratio of the observations 
and the good image quality, the differences in the 
dust emission morphology may be reasonably at- 
tributed to different radial dust properties. As- 
suming that the dust emission is optically thin (as 
discussed in § 16.21 this is a good approximation for 
most of our objects), the observed surface bright- 
ness is proportional to the emitted flux expressed 

by 

F V (R) oc H d {R) ■ k u (R) ■ T{R), (1) 

where is the dust surface density, v is the fre- 
quency of the observation, fc„ is the dust opacity at 
the frequency is, T is the dust temperature, and R 
is the orbital radius. Once corrected for the stellar 
distance, the disk inclination and the synthesized 
beam shape, different surface brightness profiles 
can result from different radial profiles of the dust 
density, opacity and/or temperature. 

On the following, we will compare the observed 
disk emission with theoretical models to (i) de- 



rive the dust properties, particularly density and 
temperature, from the available observations, (ii) 
quantify the differences between disks and (iii) in- 
vestigate the origin of the different dust emission 
morphologies. 

5. Disk model and data analysis 

To analyze the observed disk emission we com- 
pare the measured complex visibilities with a the- 
oreti cal model based on the two layer approxima- 
tion (jChiang fc Goldreich Ifl997h and on the simi- 
larity solution for the disk surface density of a thin 



keplerian viscous disk (ILv nden-Bcl 



1974t IPringle Ill98ll : lHartmann et al 



& Pringle 



199J). The 



basic properties of the model are summarized in 

SB 



The disk model that best fits the observations is 
found by adopting \ 2 as the maximum likelihood 
estimator, where \ 2 is defined by 

X 2 = E i( Re o - Rp2 t) + ( /m ' - Im 2 )] ■ w; (2) 

Re and Im are the real and imaginary part of 
the observed (suffix „) and theoretical (suffix t ) 
complex visibilities. The weight w of each data 
point in the uv plane is given by 



1 



r] s rj a Ay/2AvT a 



(3) 



where kb is Boltzmann's constant, T sys is the sys- 
tem temperature, A is the antenna area, Av is the 
band width, r acc is the integration time on source 
and rj s and rj a are the system and antenna effi- 
ciencies respectively. To minimize % 2 and evaluate 
the constraints on the model parameters we use a 
Markov Chain Monte Carlo method as described 
in Appendix [Bj 

5.1. Disk surface density 

The most common approach to derive the disk 
surface density distribution from millimeter and 
sub-millimeter observations has been to adopt a 
power law parameterization of S in the form 

(R \ p 

S( J R) = SW-^J , with R in < R < R out , 

(4) 

where Si is the surface density value at an ar- 
bitrary radius R\, and Ri n and R ou t are the 



4 



disk inner and outer radii. This parameteriza- 
tion w as initially moti vated by the empirical re- 
sults of lHavashil (1981), who derived a power law 
surface density distribution for the solar nebula 
with p = 1.5 between 0.35 and 36 AU. This sur- 
face density parameterization has been recently 
revised to accommodate a number of theoreti- 
cal and observational issues. First, there is no 
physical justification for a power law disk sur- 
face density in t erms of disk formation and evo- 



lution (see, e.g., iHueso^Gumot] [200^). More 



over, such a distribution must be artificially lim- 
ited between an inner and outer disk radius to 
obtain a finite disk mass. In addition, it has 
been demonstrated that the power law parame- 
terization fails to explain the differences in the 
radial extensions of the dust and the gas emis- 
sion that are observed in a number of interme- 



diate mass pre-main sequence stars | 


3 ietu et al. 


20051 llsella et al. 112007; Hughes et al. 


2008 


, here- 



after H08). Independently it has been suggested 
that the gas distribution in the early solar system 
would be better explained by a surface density dis- 
tribution of t he form E(i?) cx R- 1 ' 2 x exp(-i? 3 / 2 ) 



(iDavis 1 120051 



Following H08, we adopt the similarity solution 
of the surface density of a thin keplerian disk sub- 



ject to the gravity of a point mass M* (jPringle 



19811 : iLvnden-Bell fc Pringle I Il974h in the form 
presented bv lHartmann et al. I (|1998l ): 

r (2-7) 



E(r,t) 



C 



t (2-7) exp 



t 



(5) 



where C is a normalization constant, r is the stel- 
locentric distance expressed in the units of a radial 
scale factor R\ (r — R/R\), v\ is the disk viscos- 
ity at radius R\, 7 is the slope of the disk vis- 
cosity v(R) oc i? 7 , t is the non-dimensional time, 
t = t/t s + 1, t is the age of the disk and t s is the 
disk viscous time at the radius R\ defined by 



1 



1 i?? 

3(2 - 7)2 v x 

Based on Eq. [5] we demonstrate in Appendix [C] 
that i?i is the radius containing 63% of the disk 
initial mass (at t = 0). 

As pointed out by H08, this form of the surface 
density has the particular characteristic of falling 
off exponentially at large disk radii, thereby pro- 
viding sufficiently dense gas in the outermost disk 



regions to explain the observed radial extent of 
the gas emission. Moreover, since the similarity 
solution relates the surface density in Eq. [5] to the 
age of the system, we can now investigate the de- 
tails of disk evolution using our millimeter wave 
observations of the disk emission. 

The surface density as expressed by Eq. [5] in- 
cludes a significant number of unknown quantities, 
C, v\, 7, t s and R±, which cannot be constrained 
by the observations. A form more suitable for 
model fitting can be written by taking the deriva- 
tive of Eq. [5] with respect to time and introducing 
the mass flow 



M(r, t) = Ct ( 2 -7) exp 



1 - 



42-7) 



2(2- 7 )r( 2 " 7 ) 



t 



(7) 



Since disk evolution is governed by the conserva- 
tion of the total angular momentum, the disk must 
expand while matter is accreting on the central 
star. Thus M must change sign at a transition 
radius, R t , where 



R t — Ri 



t 



2(2-7) 



1/(2-7) 



(8) 



The resulting mass flow is directed inward for R < 
Rt (accretion) and outward (expansion) for R > 
Rt- 

The surface density X can be rewritten in the 
form 



x exp 



1 



2(2-7) 



(2-7) 



1 



(9) 



where we adopt the physical radius R and group 
all the other unknown quantities within S t (i.e., 
the surface density at the radius R t ). FigureH]dis- 
plays the behavior of the surface density for dif- 
ferent values of 7. When 7 = 0, Eq. [S] becomes 
a Gaussian law, while for negative 7 the surface 
density has a maximum at 



Rt. 



= R t x (-2 7 ) 1/(2 -^. 



(10) 



Our new spatially resolved observations enable us 
to constrain E t , R t and 7 and define the disk sur- 
face density distribution. 
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5.2. Disk structure and emission 

The flux emitted by the circumstcllar dust can 
be computed by solving for the structure of a pas- 
sive keplerian disk (i.e., one that is heated only by 
the stellar rad iation) adopting the two-la yer ap- 
proximation of lChiang fc Goldreich I (|lQ97t> . If the 
disk is vertically optically thick to the stellar radi- 
ation, its thermal structure is characterized by a 
surface layer temperature T s , which is appropriate 
for regions where the optical depth to the stellar 
radiation is < 1, and by a disk interior tempera- 
ture Ti, characteristic of deeper disk regions. Both 
temperatures can be calculated as function of the 
orbital radius R by iterating on the v ertical disk 



structure (see iDullemond et al. 1120011 ). Assuming 



hydrostatic equilibrium between the gas pressure 
and the stellar gravity, the disk has a flared geom- 
etry with the opening angle increasing with the 
distance from the star and the vertical gas distri- 
bution expressed by a Gaussian law normalized to 
the surface density distribution described by Eq.[9] 
Due to the radial exponential fall-off of the disk 
surface density, the very outermost disk regions 
are optically thin to the stellar radiation and the 
two-layer approximation can not be applied. For 
densities typical of TTS disks, the transition to 
this optically thin regime occurs at a radius Rd 
which is much larger than the transition radius 
Rt- Given the dust density and temperature for 
R > Rd, we adopt Rd as the disk outer radius with 
negligible effects on the strength of the overall dust 
emission. 

Once the disk thermal structure is known, the 
continuum dust emission can be computed by 
combining the flux arising from the disk interior 
Fl with the flux from the disk surface layer F*. 
These are expressed respectively by 



F'„ = 2;rcosi J \ 1 - exp 



-nR)K 



R 



xB v [Ti{R)) — dR (11) 



and 



F: 2k J \l exp 



-E(R)kt 



COS 2 



R 



(|Dullemond et al. I l200lh IChiang et al. I 120011) . 

where v is the frequency, d is the distance to 
the source, i is the disk inclination with respect to 
the plane of the sky (i = for face-on), B V (T) is 
the Planck function, AS is the column density in 
the disk surface and kl' s are the dust opacities at 
the disk mid plane and surface, as discussed in the 
next section. The disk inner radius Ri n is fixed at 
the dust evaporation distance and varies between 
~0.03 and ^0.5 AU for th e stellar luminosities 
characteristic of o ur sample ( Isella fc Natta 20051: 
llsella et al. Il2006l) . 



5.3. Dust opacity 

To calculate the disk structure and emis- 
sion we adopt the optical constants of astro- 
nomical silicates and carbonaceous materials 



x B v [T 8 (R)]AE(R)k s l/ ^dR (12) 



(|Weingartner fc Draine1l200ll : IZubko et al. Ill996f) 

The dust opacity is calculated assuming compact 
spherical grain s and adopting the fr actional abun- 
dances used bv lPollack et al. ( 1994 ) and results in 
a dust /gas ratio close to 0.01. We assume a grain 
size distribution of the form n(a) oc a~ q between 
0.01 /xm and 10 cm, where q is a free parameter 
of the model. Figure [5] shows the dust opacity at 
1.3 mm (fci.3, dashed line) and the slope f3 of the 
dust opacity (k\ oc \~@) calculated between 1 and 
7 mm (solid line), as a function of the slope of 
the grain size distribution q. ki,a reaches a maxi- 
mum value of 1.9 cm 2 per gram of dust for q — 3.9 
and decreases to values smaller than 0.4 cm 2 per 
gram of dust for q < 3.4 and q > 4.5. The slope 
(3 increases with q, and varies between [3 = 0.1 
for q = 2 (when the opacity is dominated by 10 
cm size grains) to j3 = 1.7 for q > 4.5 (when the 
opacity is dominated by sub-micron grains; see 
Natta et al. 2007 and references therein for more 
details on the variation of /3 with the grain size 
and composition). 

To solve for the disk structure, we adopt differ- 
ent values of q for the disk interior and the disk 
surface layer. To a first approximation, as long as 
the disk is optically thick to the stellar radiation 
neither the disk structure or the millimeter-wave 
dust emission depend on the assumed dust opacity 
in the disk surface layer. We therefore fix the value 
of q at 5 in the disk surface so that the opacity is 
dominated by the sub-microns grains which are 
generally required to explain the si licate features 
observed between 10 and 20 /im (jFurlan et al. 



G 



2006). The value of q in the disk interior is a 
free parameter and is chosen so as to reproduce 
the measured slope a of the spectral energy dis- 
tribution (F v cx v a ) between 0.8 and 7 mm. The 
resulting values of a, q, fci.3 and are presented 
and discussed in § 16.11 

It is important to emphasize that since we are 
analyzing the spatially resolved dust emission ob- 
served at a single wavelength, we cannot disentan- 
gle the radial variation of the dust opacity from the 
variation of dust surface density. In effect, the ob- 
servations constrain the product fci.3(i?) x E(-R) 
(see § HJ. Since the radial variation of the dust 
opacity has not been quantified observationally, 
we assume that ki.3 is constant throughout the 
disk. 

6. Results of the model fitting 

The model comparison with observations is de- 
scribed in Appendix [B] For all the objects we 
obtain good fits to the observations, with reduced 
X 2 close to 1. The dust properties adopted in the 
model fitting are presented in the Table |3J which 
list the slope q of the grain size distribution (col- 
umn 4) and the corresponding dust opacity at 1.3 
mm (column 5). Disk inclination, the position 
angle of the projected disk major axis, the tran- 
sition radius Rt, the surface density S t at Rt and 
the value of 7 corresponding to the best fit mod- 
els are shown in Table [5] The last two columns 
of Table [5] show the total disk mass and the ra- 
dius Rd at which the disk becomes optically thin 
to the stellar radiation. For each object, the ra- 
dial profiles of the disk surface density the 
cumulative disk mass Md{R), the disk surface and 
interior temperature profiles T g( j(ii) and the cu- 
mulative emission at 1.3 mm Fi^(R) are shown in 
Figure [71 The comparison between the observed 
and the best fit model visibility profiles are shown 
in Figure [H and the residuals, defined as the dif- 
ference between the observed maps of the 1.3 mm 
dust emission and the models, are shown in Fig- 
ure m 

6.1. Properties of the dust 

Figure [6] shows the spatially integrated fluxes 
of the observed disks between 0.5 and 7 mm. 
CARMA observations are represented by filled 
squares while open squares depict data from the 



literature. The solid curves show the best fit mod- 
els characterized by the values of spectral index a 
{F v cx v a ), slope of the dust opacity, slope q of 
the grain size distribution and dust opacity at 1.3 
mm fci.3 reported in column 2-5 of Table [U We 
derive values of a between 2.4 (DG Tau) to 3.5 
(LkCa 15), which lead to values of between 0.5 
and 1.7. If the disk emission is optically thin at 1.3 
mm, a and are tightly correlat ed with — a — 2 



e.g.. lBeckwith fc Sargent I ll99ll) . In practice, the 



denser regions of the disk are partially optically 
thick at 1.3 mm and > a — 2. The difference be- 
tween and a — 2 depends on the ratio of optically 
thick to optically thin emission from the disk and 
is only few percent in CY Tau, DN Tau, UZ Tau E 
and GSS39. However, it increases to about 20% 
in SR 24 S and DG Tau, where the disks are opti- 
cally thick to the 1.3 mm emission within a radius 
of 17 AU and 25 AU respectively. 

The dust opacity fci.3 obtained by the SED fit- 
ting strongly depends on the assumed grain size 
distribution discussed in § 15.31 Actually, we can 
reproduce the observed spectral indexes a with 
very different dust opacities if we fix the slope of 
the grain size distribution q and keep a max as a 
free parameter. Values of a max and fci.3 for the 
case q = 3 are shown in the last two columns of 
Table |U This latter grain size distribution leads to 
dust opacities larger than the case with fixed a max 
and variable q by a factor 2-20. As a consequence, 
much smaller surface densities and disk masses are 
required to reproduce the observed dust emission. 

Independently on the assumed grain size distri- 
bution, values of < 1 imply that the dust opacity 
is strongly influenced by dus t grains larger than 
1 mm (e.g. iNatta et al. 1120071) . This suggests that 
the circumstellar dust around CY Tau, DG Tau, 
DN Tau, DR Tau, RY Tau, UZ Tau E, GSS 39 and 
TW Hya has undergone important grain growth 
processes. However, the small value of the spec- 
tral index for DG Tau can now be probably ex- 
plained by the strong contribution of the optically 
thick emission at 1.3 mm. By contrast the dust 
properties in LkCa 15 and GO Tau are more simi- 
lar to that found in the ISM. Moreover, disks with 
small value of tend to be characterized by small 
value of 7 and therefore to steeper surface density 
profiles for R > R t . This suggest that the dust 
surface density may be correlated with the dust 
properties and in particular that small and large 
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dust grains might have different radial density dis- 
tributions as suggested by recent theoretical mod- 
els of growth and radial migration of dust g rains 



in proto-planetary disks (jBrauer et al. 1 12008). 



6.2. Radial profile of the surface density 

In the panels of the first and second columns 
of Figure [7] we show for each object the surface 
density and the cumulative flux at 1.3 mm (solid 
line), the 1-cr uncertainty range (shaded region) 
and the spatial resolution provided by our inter- 
ferometric observations (vertical dashed line) . The 
dust emission coming from spatially resolved disk 
regions varies from about 85% (CY Tau, LkCa 15, 
MWC 275, GSS 39) to less than 30% (DR Tau, 
SR 24 S) of the total flux. As a consequence, the 
constraints on R t , 7 and S t vary from few percents 
to about 30-40%. 

The derived transition radii range from 17.5 
and 110 AU. For the assumed dust opacities and 
stellar distances, the radius Rd at which the disk 
becomes optically thin to the stellar radiation - in 
effect the disk outer radius (§ 15. 2|) - varies from 
73 AU in the case of TW Hya to 670 AU for 
GO Tau. These values are close to the disk outer 
radii inferred from the analysis of th e optically 
thick CO emission ( Simon et al. 2000h . but larger 
than existing determinations of the disk outer 
radii based on power law surface density models 



(|Andrews et al. 1120071 ). In the case of MWC 275, 



we find Rd — 520 AU. This radius is in very good 
agreement with the gas extent inferred from the 
CO line emission but it is twice as large as the disk 
outer radius implied by the same observations if a 
power law surface density disk model is adopted 
( Isella et al. ~l l2007l) . t appears that there is in fact 
no discrepancy between the radial extents of the 
gas and dust emission. The outer radii determina- 
tion of the dust extent based on the exponential 
fall-off for the surface density are quite compara- 
ble with the radii derived from observations of the 
CO emission (see also H08). 

From Eq. [9l it is clear that the disk surface den- 
sity is characterized by 7, which ranges from -0.8 
to 0.8. For R> R t , radial density profiles become 
steeper with decreasing 7, as shown in Figure [5] 
This is further illustrated by the plots of £ in the 
panels of the first column of Figure [71 In Fig- 
ure EH LkCa 15, CY Tau, DC Tau and TW Hya, 
have negative values of 7 with a probability of at 



least 68% (1-cr, see Appendix [Q. From Eq. [10] 
(§ I5.1|) . the corresponding surface densities have 
maxima at 89, 41, 21 and 13 AU respectively. For 
LkCa 15, the CARMA observations clearly resolve 
the surface density maximum and the almost flat 
profile of £ between 50 and 100 AU. For the other 
objects the predicted maxima lie inside the reso- 
lution of our observations. Disks with values of 
7^0 are shown in Figure I7b| and for 7 > in 
ITcl For almost all disks, £ is well constrained for 
R > 40 AU. For DR Tau and SR 24 S, almost 80% 
of the observed emission arises from the innermost 
spatially unresolved disk region so that the disk 
structure is poorly constrained. The structure of 
GO Tau disk is also poorly constrained due to the 
low S/N ratio of the observations. 

The inferred disk structure depends weakly on 
the adopted grain size distribution if the SED is 
used to constrain the dust opacity. In fact, the disk 
mid plane temperature Ti(R) varies by less than 
5% between the two different grain size distribu- 
tion models shown in Table |4j Consequently, since 
£(-R) oc T _1 (i?), 5% is also the maximum varia- 
tion observed in the profile of the surface density. 
We obtained model fits for a subset of the disks 
using the alternative grain size distribution listed 
in Table [H and have verified that the variations 
on 7 and R t are much smaller that the respective 
uncertainties. 

6.3. Disk mass 

The radial integration of the disk surface den- 
sity leads to the cumulative disk mass M(R) pre- 
sented in the panels in column 3 of Figure [7J and 
to the total disk mass Md in Table \5[ For the 
adopted grain size distribution and dust / gas ratio 
(§ 15. 3p . Md varies by more than a order of magni- 
tude from ~0.02 to ~0.4 M Q , with uncertainties 
between 20% and 200% depending on how well £ 
is constrained. Md is. In the case of DG Tau, 
the mass of the disk is comparable with the stellar 
mass. Note however that stellar masses derived 
using DM97 stellar evolution models are probably 
a lower estimate of the real stellar mass. Indeed, 
the disk mass is onl y 35% of the stellar mass if 
Baraffe et al. I (|l998h models are used to derive 
the stellar mass from the HR diagram (see Ap- 
pendix [X]). Disk mass is 15% of the mass of the 
parent star for DN Tau and GM Aur, while it is 
only a few percent for all the other objects. 
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Disk masses strongly depend on the dust opac- 
ity at 1.3 mm which may vary by a large factor if 
a different grain size distribution is adopted. The 
case q = 3.0 discussed is Section loTTl leads to disk 
masses between a factor 2 and 20 smaller than 
what discussed above. This introduces an addi- 
tional large uncertainty on the disk masses derived 
from millimeter-wave observations. 

As shown in the cumulative mass plots, most of 
the disk mass is concentrated in the outermost disk 
regions, independent of the exact surface density 
profile. By contrast, since the disk scale height 
H(R) increases with distance from the star, the 
dust+gas volume density in the disk mid-plane 
[po oc Ti(R)/H(R)] is at a maximum close to the 
inner disk radius. A comparison of Ma with the 
minimum mass of the solar nebula (0.02 Mq, see 
Figure [7]) shows that this amount of material is 
contained within radii ranging from — 10 AU (for 
DG Tau) to -60 AU (for DM Tau, MWC 275). 
These values are not too far from the minimum 
mass solar nebula outer radius of 36 AU postu- 
lated by Hayashi (1981). We will return to this 
point in § 17.31 

6.4. Disk temperature 

The panels in the last column of Figure [7] dis- 
play the radial profile of the dust temperature in 
the disk interior, Ti (thick line), and in disk sur- 
face layer, T s (thin line). As discussed in § 15.21 T s 
scales roughly as i? -1 / 2 due to the dilution of the 
stellar radiation and to the variation of the Planck 
mean dust opacity. It varies between —1500 K at 
the disk inner radius, and 10 to 40 K at Rd, where 
the disk becomes optically thin to the stellar ra- 
diation in the vertical direction. Since the dust 
column density in the disk atmosphere is a tiny 
fraction (— 10~ 4 ) of the total disk surface density, 
the contribution to the observed millimeter emis- 
sion from the disk atmosphere is negligible (see 
Eq.rJU). 

Inside a radius of 50-100 AU, the dust tem- 
perature Ti in the deeper disk region also scales 
as i? -1 / 2 , as expected for an irradiate d optically 
thick disk (jKenvon fe Hartmann 1119951 ) . At larger 
distances the dust density is low enough that the 
disk interior becomes progressively optically thin 
to the emission of the disk atmosphere (which 
heats up the disk interior) and to its own ther- 
mal emission (which cools the disk interior). The 



resulting effect is that the temperature profile de- 
viates from the i? _1//2 relation. The deviation is 
larger in disks characterized by small values of 7 
and a correspondly steeper decrease in the surface 
density (see Figure ITal and [7b)). In these objects, 
the temperature of the disk region that dominates 
the observed millimeter emission assumes almost 
constant values between 20 and 40 K. By contrast, 
Ti decreases monotonically in disks with the high- 
est 7 (i.e., DM Tau, GO Tau and UZ Tau E). In 
no cases does the disk temperature fall below 10 K 
which is generally assumed to be the equilibrium 
temperature with the interstellar radiation field. 

6.5. Disk orientation 

Disk inclinations and the position angles de- 
rived from the observations are uncorrelated 
and randomly distributed as expected in ab- 
sence of a preferential disk orientation in space. 
For CY Tau, DG Tau, DN Tau, DM Tau, 
LkCal5, UZ Tau E, GM Aur, GSS39, SR24, 
TW Hya and MWC 275 both inclinations 
and position angles are in agreement wit h pub- 



lished values within 2a ( Jensen et al. _ 1996; 



Dutrev et al. 



2004; 



llsella et al 



Pietu et al. 
20071 ) 



1997E ISimon et al 
20061: 



2000 



Qi et al. 



Andrews et al. 



20071: 



For GO Tau, DR Tau and 



RY Tau we derive disk inclinations of 25°±25°, 
37°±3°and 60°±3°respectively, consid erably lower 
than t he 66°, 72° and 86° suggested bv lAndrews et aL 
(|2007l ). On the other hand our RY Tau re- 
sults agree well with recen t optical observations 



(jSt-Onge fe Bastien I [2008) suggesting that the 



di screpancies may be du e to the lower S/N ratio 



of I Andrews et al. I (|2007l ) observations compared 



with the new CARMA results. 

7. Discussion 

From the results of § we conclude that the di- 
versity in the surface brightness profiles observed 
in our 1.3 mm dust continuum maps (§ [4j is due 
mainly to different surface densities, S(i?), with 
the disk temperature Ti(R) having only a minor 
effect. Of course, this assumes that the dust opac- 
ity is constant with radius as discussed in § 15.31 
Most of the discussion below is therefore devoted 
to examining the dust density profiles and inves- 
tigating the possible sources of the observed vari- 
ations. 



9 



7.1. Disk evolution 

Figure \W\ compares the disk transition radius 
Rt and the stellar age (from Appendix [5J ■ Over 
the ~5 Myr span of our sample ages, Rt appears 
to increase with stellar age from about 20 to 100 
AU. Applying the non-parametric Spearman rank- 
order test (see, e.g., Press et al. 2007), the cor- 
relation coefficient is r = 0.42 and the proba- 
bility that the data are randomly distributed is 
12%. However, limiting the analysis to the ob- 
jects in Taurus-Auriga (filled squares) the corre- 
lation coefficient increases to 0.98 and the proba- 
bility of a random distribution falls below 0.1%. 
If R t has a power law dependence of the form 
R t — Ro + C -t v , where t is the stellar age in Myr, 
we find that R t increases as \fi , with r\ = 0.5±0.4, 
Ro = 17 ± 10 AU and C = 37 ± 20. As shown 
in Figure 1111 this disk expansion is accompanied 
by a decrease in the mass accretion rate roughly 
described by M acc oc f~ 1>4±0 " J . In this case, the 
Spearman test indicates that the probability that 
the data are randomly distributed is 2% (r=-0.62) 
both for the full sample and for the objects in 
Taurus- Auriga. As illustrated in Fig. [T^l there 
appears to be no correlation between 7 and stellar 
age (r=0.18). 

The observed increase of the disk transition ra- 
dius, R t , and decrease of the mass accretion rate, 
M a cc, suggests that the our sample disks repre- 
sent an evolutionary sequence. Stars with different 
mass, spectral type and luminosity are accompa- 
nied by disks whose characteristics seem to vary 
similarly with time, i.e. they evolve in a similar 
way. Can we interpret this in terms of the viscous 
disk model of § 15.11 .'' 

To answer this question we can start from Eq. [8] 
and express the temporal variations of Rt as: 

1/(2-7) 



Rt — R-\ 



1 



t 

- + 1 



(13) 



.2(2-7) U 
and from Appendix |Pl the mass accretion rate 

M„(0) ft x-(./ a -7)/C-T) 

(14) 

where 7 is assumed to be constant with time. The 
values of 7 derived for our sample (Table [5]) lead to 
R t oc t - 3 - - 8 and M acc cx H 1 - 2 " 1 - 4 ). The agree- 
ment with the relation derived above from our ob- 
servations suggests that the values of R t and M acc 



result from viscous evolution of disks formed with 
similar masses and radii over a time interval of 
about 5 Myr. 

7.1.1. Initial disk properties and time scale for 
the disk evolution 

The initial disk radius, R±, the initial disk mass 
Md(0) and the viscous time scale at R\, t s , can be 
estimated from Eqs. [T3land[14l Assuming values 
of 7 between -0.8 to 0.8, we obtain R x =25-40 AU, 
M d (0) = 0.05-0.4 Mq, and t s =0.1-0.3 Myr. In 
Appendix 0we showed that R\ is the radius con- 
taining 63% of the initial disk mass M d (0), while 
about 90% of the initial disk mass M d (0) is con- 
tained within about 2R\ — 50-80 AU. As noted in 
Appendix [3 stellar ages increase by a factor 3-10 
if B98 models are used instead of DM97. Never- 
theless, the transitional radii, R t , and the mass 
accretion rates, M acc , still correlate with the stel- 
lar ages with similar correlation coefficients. R\ 
and Mrf(0) are influenced very little by different 
evolutionary models but the resulting viscous time 
scales are about a factor 4 larger. In the following 
discussion we focus on deriving constraints on the 
processes that govern the disk viscosity and the 
disk formation from R\, M^(0) and t s . 

7.1.2. Implications for disk viscosity 

From the disk radius R\ , the corresponding vis- 
cous time scale t s and the parameter 7, we can de- 
rive the absolute value of the disk viscosity and its 
radial profile v(R). From § 15. R v(R) is expressed 

by 

I R\ ( R 



u(R) 



(15) 



3(2-7) 2 t s 

Since different values of 7 are derived from our ob- 
servations of the millimeter dust emission, the disk 
viscosity is probably characterized by a variety of 
radial profiles. Thus for 7 > (DM Tau, GO Tau, 
UZ Tau E, GM Aur, GSS 39 and MWC275) the 
viscosity increases with radius. For 7^0 (CY 
Tau, DN Tau, DR Tau, RY Tau, SR 24) the vis- 
cosity is virtually constant with radius, and for the 
few objects with 7 < 0, it decreases. 

The viscosity in a disk is generally attributed to 
some degree of turbulence that may originate from 
different physical processes. Two possible sources 
of tur bulence are the magnet o-rotational insta- 
bility (jBalbus fc Hawlev Ill991r ) and gravitational 
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instability (e.g.. iLodato fe Rice 112004 ). Magneto- 
rotational instability (MM) requires weak mag- 
netized disks and a minimu m ionization frac- 
tion of about 10" 13 at 1 AU (jBalbus fc Hawlev 
2000). For the disk to be MRI active the max- 



imum surface density must lie between 10 and 
100 g/cm 2 , depending on the ionization source 



(jTerquem et al. 1 12008). This condition is satis- 
fied in most of our objects (see Figure [7|). For 
DG Tau, SR 24 S and RY Tau, however, the 
high density at the disk mid-plane inside a ra- 
dius of about 30 AU probably prevents ionization. 
In such regions, where the MRI cannot operate, 
viscosity may originate from gravitational insta- 
bilities (GI). Classically GI can be parameterized 
by the Q-value (Toomre 1964) 



Q = 



7rGS(i?) 



(16) 



where 51 is the keplerian angular velocity (f2 = 
■J GM+/R 3 ) and c s the sound of speed. When 
Q < 1.5 the disk is gravitational unstable and de- 
velops spiral waves to transport angular momen- 



tum outward and mass inward (Lodato fc Rice 



2004). Adopting Eq. IE2I for the sound speed in 
the disk, we can rewrite Q in the form 



Q ~ 230 



Q.5M G 



1/2 



R 



10AU 



-3/2 



E(ii)- 1 r i (i?) 

(17) 

where as usual £ and Ti are the disk surface 
density and interior temperature as in Figure [7] 
Among the observed objects, DG Tau has a grav- 
itationally unstable disk between ~ 20 and ~ 60 
AU (solid curve in Fig. [T3]). For LkCa 15 (long 
dashed curve), Q is close to 1.5 between 60 and 
160 AU while in all the other objects (short dashed 
curves) Q is well above the instability threshold. 
Since Q oc E(i?) _1 , it strongly depends on the as- 
sumed grain size distribution as discussed in § 16.11 
In the case of the grain size model with q = 3 
and variable a max , Q is always above the 1.5 and 
gravitational instabilities cannot occur. However, 
the smaller surface density leads to an higher ion- 
ization fraction, facilitating the role of the MRI. 
It seems therefore that both MRI or GI may be 
applicable in the observed objects. 

We compare the viscosity profile v(R) from 
Eq. [15] with the theoretical expectation of MRI 
and GI models by adopting the classical a viscous 



1/2 



disk parameterization (Shakura & Sunyaev 1973) 
and writing the stress parameter a in the form (see 
Appendix [E} 



a(R) 



3 x 10~ 2 / t s 



(2- 7 ) 2 \lMyr 



R 



10AU 



(7-3/2) 



0.5M 



7 r, y 2 -^ 

VIOAuJ 

1/2 

T-\R). (18) 



It is of some importance to note that a is generally 
assumed to be constant in analytic modeling, al- 
though this choice has no physical justification and 
numerical simul ations show variations of a b oth in 
space and time ( Nelson fc Papaloizou 1 12003). For 
a simple disk model that assumes Ti(R) cx i? -1 / 2 
(see §[6]), constant a corresponds to the case 7 = 1. 
However, since the disk temperature T deviates 
from the R- 1 ' 2 profile (§[674]) and we observe cases 
of 7 < 1, a cannot be constant in the observed 
disks. Figure [T4l shows the behavior of a for 7 > 
(upper panel), 7 ~ (middle panel) and 7 < 
(lower panel). In all cases, a decreases with the or- 
bital radius R and it may vary by almost 4 orders 
of magnitude between 0.5 and 10~ 4 . For 7 > 0, 
a has a shallow dependence on radius and ranges 
from about 0.03, at 1 AU, to about 0.005, at 100 
AU. For 7 < 0, a varies more rapidly and assumes 
values larger than 0.1 inside radii of few AU and 
smaller than 0.001 outside about 30 AU. 

Numerical simulations of keplerian disk af- 
fected by MRI suggest that a can ran ge from 



Balbus 



0.005 to 0.6 (see, e.g., the review of 
20031) . Although detailed comparison of the MRI 
model and our observations is beyond the scope of 
this work, we note that a decrement of a with 
the radius can be achieved for particular val- 
ues of the magnetic field st rength and geometry 
(jPapaloizou fc Nelson 1120031 ). As discussed above, 
the DG Tau disk can be gravitationally unstable 
between 20 and 60 AU. In this region we calcu- 
late values of a between 0.001 and 0.003 which 
are slightly smaller than a ~ 0.05 predicted by 
numer ical simulation of gr avitationally unstable 
disks (|Lodato fc Rice 112004) ■ 

Although our analysis is very qualitative, MRI 
seems to be able to account for the values of a 
derived in our sample disks, and, perhaps even 
the radial profiles. Nevertheless, it is still un- 
clear whether circumstellar disks can be com- 
pletely MRI active, or if the MRI is effective 
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only in the very inner part of the disk (e.g., 
Chiang fe Murrav-Clav ll2Q07h . More importantly, 
it is still an open question why MM would oper- 
ate in different ways in our sample disks, leading 
to the variety of values of 7 and a discussed above. 

7.1.3. Implications for parent cores 

Our viscous time scales and stellar ages sug- 
gest that the assumption t t s made in § 15.11 is 
probably not appropriate for the youngest objects, 
e.g., DG Tau, DR Tau, GSS 39 and SR 24 S. For 
these objects, with little time for disk expansion, 
the dust radial distribution should trace the ini- 
tial disk structure resulting from the collapse of 
the parent core and provide insight into as yet un- 
known details of the the formation process (e.g., 
Andre et alll2008l ). 

A simple relation between circumstellar disks 
and core properties can be obtained on the as- 
sumption that disks for m from the collapse of 
rigidly rotating cores ([Goodman et al. 1993; 
Shu et al. Hl977l L If w c is the core angular velocity 
and if the core angular momentum J is conserved 
during the collapse, the disk initial radius R\ can 
be expressed as: 



core density gr adient of p(r) oc r 2 and a core ra- 



i?i ~ 25 



(lO- 14 s-0 ( 



AU 



(19) 



(|Hueso fc Guillot 1120051: iDullemond et al. I I2OO6). 



Here we have assumed that (i) the disk centrifu- 
gal radius R c (i.e, the radius at which the angular 
momentum of the disk is equal to the angular mo- 
mentum of the parent core) is similar to the radius 
that contains about 90% of the initial disk mass 
(i.e, R c ~ 2Ri, Appendix [Cj, (ii) the disk mass is 
negligible compared to the mass of the central star, 
(iii) the disk does not significantly expand while it 
is still accreting material from the e nvelope, (iv 



the tempera ture of the core is 10 K (jKirk 112007 



Jiiina 1999f) and (v) the magnetic field does not 
play a significant role in the core collapse. 

For our stellar masses (Tabled]) and disk initial 
radii between 25 and 40 AU f§ l7.1|) , we derive core 
angular velocities between 5 x 10~ 15 and 2 x 10 -14 
s _1 . Assuming simple radial profiles of the core 
density, we can estimate the specific core angu- 
lar momentum j (i.e., the angular momentum per 
mass unit j = J/M) required to form disks with 
initial radii in the observed range. We assume a 



dius of 0.05 pc (jCaselli et al. 1120021 and references 
therein). The core specific angular momentum, 
j = (2/9)w c i? 2 ore , then lies between 8 x 10~ 5 and 
4 x 10~ 4 km/s pc. 

For dense core, specific angular momenta j = 
(0.5 — 4) x 10~ 3 km/s pc have been derived from 
the measurements of velocity gradients of chemical 
tracer such as NH3 and N2H+, and for cores with 



mass e s in the range 1 to 1 Mq ([Goodman et al. 



1993; ICaselli et al. 1 120021 ). These values are an 



order of magnitude larger than required to form 
disks with initial radii between 25 and 40 AU, sug- 
gesting that about 10% of the specific core angular 
momentum and about 30% of the core mass are 
conserved during the disk formation. 

7.2. Interpreting "transitional" disks. 

The surface density in disks with 7 < (i.e., 
LkCal5, TW Hya and DG Tau) increases with the 
orbital radius and reaches a maximum at R ~ Rt 
(§[6]). This effect can mimic the inner disk clearing 
advocated to explain the deficit in the near and 
mid-infrared excess over the stellar photosphere 
observed in "transitional" disks ( Strom et al 



1989) where, it is postula ted, planetary system 
formation may have begun (jEspaillat et al. 2008). 
Thus, for LkCa 15 and TW Hya, we predict a 
gradual decrease of the surface density inside radii 
of ~60 and ^17 AU respectively, in qualitative 
agreement with the radii of the dust depleted disk 
region already inferred from mm-wave obs erva- 
tions (|Pietu et al. 1120061 : iHughes et al. 1120071 ). 

Figure [E] shows the observed spectral energy 
distribution of LkCa 15 (points) and the disk 
model of Table [5] that fits our 1.3 mm observations 
(solid line). The disk surface density is character- 
ized by 7 = —0.8 and extends up to the disk inner 
radius without any discontinuity (see Figure I7a,|) . 
At about R = 0.1 AU, the dust sublimates, form - 



ing a "puffed-up" inner rim (jlsella fc Natta 1 12005) 



which emits mostly in the near-infrared between 1 
and 3 /zm (long-short dashed lines). For R > 0.1 
AU, the optically thin disk surface layer emits in 
the mid and far infrared (short-dashed line) while 
the colder disk midplane dominates the emission 
at longer wavelengths (long-dashed line). The 
model reproduces well the observed SED. 

Until now, photoevaporation, the presence 
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of a planet, grain growth and inside-out MRI 
clearing have been the main processes invoked 
to explain the properties of transitional disks 



(120051) differs from ear l ier studies dKusaka 
Weidenschilling 1977 : Havashi 198lh and sug- 



dAlexander et al. l2006tlCalvet ll200llEspaillat et al . 
2008tlDullemond k. Dominik l2005HChiang k. Murrav-Cla"it iscribed by 



gests that the surface density distribution in the 
solar nebula at 10 5 -10 6 yr from its formation is 



20071 ). Here we suggest disk viscosity as an addi- 
tional means of producing a surface density pro- 
file and spectral energy distribution consistent 
with this class of disks. In fact, the similarity 
solution for the disk surface density predicts par- 
tially depleted inner disks whenever the viscosity 
v{R) oc R 1 decreases with the orbital radius. 

If transitional disks are explained by the sim- 
ilarity solution for the disk surface density, the 
surface densities at R < R t and R > R t are then 
tightly correlated by Eq. [9] In particular, since 7 
must be negative, the surface density must fall-off 
very quickly for R > R t (see Figure SJ. The outer 
disk in transitional objects must be therefore char- 
acterized by a rapid decrease of the surface bright- 
ness measured at millimeter wavelengths. In our 
sample this hold true for LkCal5 and TW Hya but 
it is not verified for DM Tau and GM Aur. For 
the latter two objects, the observed dust contin- 
uum emission at R > 45 AU leads to positive value 
of 7, implying that the surface density in DM Tau 
and GM Aur increases smoothly up to the inner 
disk radius (see Figure [7}. This contrasts with the 
observations of dust depleted inn er disks within 
about 7 and 20 AU (jCalvet 1120051 ). 

We suggest that this interpretation of "transi- 
tional disks" is applicable when the dust-depleted 
inner region occupies a significant fraction of the 
disk extent. In LkCal5 and TW Hya, R t is in 
fact more than 20% of the disk radius Rd (Ta- 
ble [5]). By contrast, when the dust-depleted inner 
region is only a few percent of the disk radius, 
as for DM Tau and GM Aur, the presence of a 
planet, the inside-out MRI clearing or the photo- 
evaporation by the central star are more probable 
explanations for the inner disk clearing. 

7.3. Similarity with the solar nebula 

A fundamental question regarding pre-main se- 
quence disks is whether they will evolve into plan- 
etary systems similar to our own. Addressing 
this question is difficult mainly because we have 
a limited knowledge on the properties of the so- 
lar nebula. A recent re-analysis of the distribu- 
tion of solid bodies in the Solar system by iDavis 



£(i?) ~ 1.14 x lO^-^e- 024 * 3 



(20) 



The numerical constants are chosen to re cover the 
total mass of the solar nebula, 0.02 M (IKusaka 



1973), and its specifi c angular m omentum j 



3.7 x 10" 6 km/sec pc dCox Il2000h . 

The solar nebula surface density of Eq. [20] is 
based on the assumption that 7 = 0.5, and cor- 
responds to the similarity solution expressed by 
Eq. |5] for R t = 6 AU and E t = 340 g/cm 2 . The 
total mass for the solar nebula are quite similar to 
the values measured for GSS 39, GM Aur, DM Tau 
and GO Tau. However, the solar nebula transi- 
tion radius, Rt, is at least a factor of 4 smaller 
than what found for objects with ages of 10 5 -10 6 
yr (see § IO and § 17X31) . 

Indeed, among our disk sample, only TW Hya 
may be a good match to the solar nebula, being 
the only old object characterized by a small tran- 
sition radius. If TW Hya underwent viscous evolu- 
tion similar to that discussed in § 17.11 it probably 
formed with an initial transition radius of 4-8 AU, 
similar to that predicted for the solar nebula. 

8. Conclusions 

We presented high angular resolution (0.7") in- 
terferometric observations of the 1.3 mm contin- 
uum emission from 14 pre-main sequence circum- 
stellar disks. The disk surface brightness is char- 
acterized by a range of radial profiles. Adopt- 
i ng the similarity soluti on for disk surface density 



(jHartmann et al . 1119981) and a self consistent disk 



emission model ( Isella et al. 112007 ). we derived for 
each disk a surface density radial profile defined 
by the transition radius Rt , the surface density S t 
and the slope of the disk viscosity 7, as well as 
a radial temperature profile and its orientation in 
space. Assuming a constant dust opacity through- 
out the disk, we find that the different surface 
brightness profiles are mainly due to differences 
in R t and 7, which in turn imply different disk 
surface densities. 

From a comparison of the disk surface density 
and the stellar properties, it appears that the disk 
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transition radius Rt is correlated with the stellar 
age and increases from ^20 to ~ 100 AU over 
about 3-5 Myr. This disk expansion appears to be 
accompanied by a decrease in the mass accretion 
rate. We argue that these temporal variations of 
the disk radius and the mass accretion rate sup- 
port a scenario in which disks form an evolutionary 
sequence. The observed evolution is in qualitative 
agreement with that of viscous disk models if the 
initial disk masses are between 0.05 and 0.4 M Q 
and the initial disk radii range from 20 to 45 AU. 
Note however that disk masses can vary by one or- 
der of magnitude depending on the assumed dust 
opacity. The temporal variation of the disk radius 
and mass constrains the viscous time scale t s at 
the disk initial radius R\ to be between 0.1 to 0.3 
Myr. 

The viscous disk model assumes that the disk 
viscosity v(R) scales as R 1 . Among our sample, 
7 ranges from -0.8 to 0.8, leading to a large vari- 
ety of viscosity radial profiles. Parameterizing the 
disk viscosity in terms of the stress parameter a, 
we show that a scales with radius R roughly as 
i? 7_1 . Since 7 is always smaller than 1, a must 
decrease with the orbital radius. We suggest that 
a may vary by almost 4 orders of magnitude be- 
tween 0.5 and 10 -4 . These values are in general 
agreement with MRI models, which is probably 
the main source of viscosity at the surface densi- 
ties of our disks. However, its still an open ques- 
tion why MRI should operate differently among 
our sample, leading to a large range of values of 7. 

The ages of the younger stars in our sample 
are comparable to the viscous time scale of 0.1-0.3 
Myr. It seems likely that, with little time for the 
disk expansion, the dust radial distribution should 
trace the initial disk structure resulting from the 
collapse of the parent core. For typical assump- 
tions on the core rotation, radial density profile 
and radius, we derive the core specific angular mo- 
mentum j — J/M required to form disks with ini- 
tial radii between 25 to 40 AU. We argue that j 
must range from 8 x 10~ 5 to 4 x 10~ 4 km/s pc, 
and suggest that this corresponds to about 10% of 
the specific angular momentum measured in dense 
cores. It seems therefore likely that during the star 
formation process, about 10% of the core angular 
momentum is transferred to the circumstellar disk. 
Alternatively, 10% efficiency in conserving angu- 
lar momentum implies that about 30% of the core 



mass is used to form the central star. We believe 
that the attempt to correlate the properties of cir- 
cumstellar disks and dense cores, though still very 
qualitative, is an important step forward to un- 
derstand the role of disks in conserving the angu- 
lar momentum during the star formation process. 
Clearly, the capability to investigate this aspect in 
detail are actually hampered by the sensitivity and 
resolution of the existing interferometers, which 
limit the analysis to a few bright objects. How- 
ever, the ongoing improvement of facilities such 
as CARMA and VLA, and the advent of ALMA, 
will enable us to expand the number of spatially 
resolved disks by a large factor and have a more 
complete view on the relations between dense core 
and young disks. 

We point out that in disks with 7 < 0, the 
surface density increases with radius R and 

reach a maximum at about the transitional ra- 
dius Rt- We argue that this particular behavior 
of can mimic the inner disk clearing ad- 

vocated to explain the dominant characteristic of 
some "transitional disks" , namely a deficit in the 
near and mid-infrared excess over the stellar pho- 
tosphere and the presence of an "hole" in the 
surface brightness observed at millimeter wave- 
lengths. For LkCal5, we show that a surface den- 
sity with 7 = —0.8 that extends without any dis- 
continuity up to the disk inner radius located at 
few stellar radii can reproduce both the SED and 
the 1.3 mm continuum emission. By contrast, we 
find no clear explanation in terms of the similar- 
ity solution of the surface density for the dust de- 
pleted inner disks around DM Tau and GM Aur. 
It seems likely that "transitional disks" may orig- 
inate from a large variety of effects. 

Finally, it appears that most of the disks in 
our sample are very different from the currently- 
accepted view of the solar nebula. While most of 
them have masses similar to the minimum mass 
solar nebula, 0.02 M©, their transitional radii are 
at least a factor 4 larger than the value Rt ~ 6 
AU derived from the actu al distributio n of solid 
bodies in the Solar system ( Davis 2005^ 1 . The ex- 
ception is TW Hya which has a very small disk 
radius, R t ~ 17 AU, compared to its large age, 
t = 7 Myr. We argue that TW Hya disks proba- 
bly formed with R t in the range from 4 to 8 AU, 
and may well reflect the properties of our early 
solar system. 



14 



We are indebted to Meredith Hughes and David 
Wilner for providing the SMA data of GM Aur 
and TW Hya. We thank Antonella Natta, Sean 
Andrews and Leonardo Testi for the useful dis- 
cussions and the referee for the very useful com- 
ments. We thank the OVRO/CARMA staff and 
the CARMA observers for their assistance in ob- 
taining the data. We acknowledge support from 
the Owens Valley Radio Observatory, which is 
supported by the National Science Foundation 
through grant AST 05-40399. This work was per- 
formed in part under contract with the Jet Propul- 
sion Laboratory (JPL) funded by NASA through 
the Michelson Fellowship Program. JPL is man- 
aged for NASA by the California Institute of Tech- 
nology. 



15 



A. Stellar ages and masses 



Figure [TB"| shows the l ocatio n of our sample stars in the H-R di a gram with the evolutionary models of 



D'Antona fc Mazzitellil <|l997l hereafter DM97) and iBaraffe et al. I (|l998l B98). The uncertainties on the 
stellar temperature correspond to half a spectral type, while the stellar luminosity is uncertain by 30%. 
From DM97 models, the stellar ages range from 0.1 Myr (DR Tau and DG Tau) to 7 Myr (TW Hya). 
In Taurus-Auriga alone the age spread is 1-3 Myr. However, as illustrated in the upper panel of Figure 
[TTl these ages are smaller than those inferred from B98 by factor of 2.6-10. Similarly, the lower panel of 
Figure [T7] shows that the stellar masses inferred from B98 are systematically large r than those from D M97 
by factors of 1.2-3.4. Neither model reproduces the dynamical masses derived by ISimon et al. I (|2000h for 
DM Tau, GM Aur, CY Tau and LkCa 15. Typically the DM97-derived masses are lower than the dynamical 
masses by 10-40% while B98-masses are lower by 20-30%. For the purpose of our analysis we have adopted 
DM97 models because they enable us to derive masses and ages for all the stars in our sample, while RY Tau 
and MWC 275 are outside the temperature and luminosity range of the BH98 model. 



B. Fitting process 

For fixed stellar parameters, dust opacity and disk inner radius, the dust emission model is defined by the 
state x{i, PA, R t , St, 7} where i is the disk inclination, PA is the disk position angle measured from North 
to East, Rt is the disk transition radius, S t is the disk surface density at R t and 7 defines the shape of the 
disk surface density as discussed in § 15.11 The model that best fits the observations is found through the 
minimization of \ 2 (Eq. EJ. We adopt a Bayesian approach in which the joint probabilities for the observed 
data d and the given model state x are described as the product of the probability of the observed data d 
given the model parameters x (i.e., the likelihood), and a known prior probability distribution function p(x) 
of the model parameters: 

p(x\d) oc p(x)p(d\x). (Bl) 

In this framework, the best-fit model corresponds to the state x that maximizes the a posteriori distribution 
p(x\d). 

To characterize the posteriori distribution we adopt a Markov Chain Monte Carlo (MCMC) method (see, 
e.g., Ford 2005, Fitzgerald et al. 2007). This method generates a chain of states x sampled from a desired 
probability function p(x), whose equilibrium distribution is equal to the posteriori distribution p(x\d). In 
our specific case, since complex visibilities errors are described by a normal distribution (Wrobel & Walker 
1999), the probability of the observables (i.e, the real and imaginary part of the complex visibility) given 
the model state x, is expressed by the x 2 distribution. To a first approximation this is proportional to 
e -x (a;)/2_ j n the most general case, where all model states have the same prior probability and p(x) ~ 1, 
the posteriori joint probability is also roughly proportional to e~ x ( x )/ 2 . The state that maximizes p(x\d) 
therefore correspond to the state that minimizes x 2 (cc), justifying the choice of the \ 2 as the maximum 
likelihood estimator. 

The MCMC fitting is realized following iFord I (|2005l ) and is briefly summarized here. First, the chain 



is constructed using the Metropolis-Hastings algorithm with the Gibbs sampler, and generating a chain of 
models where the (n+l)-th state depends only on the n-th state through a specified transition probability 
g(n|n+l). In practice, at each n-th step of the chain we (1) generate a random trail state x 1 adopting 
a transition probability described by a Gaussian distribution centered on x, (2) calculate x 2 (a;'), and (3) 
accept the trail state as the new (n+l)-th state of the chain if it satisfies an acceptance probability a(x'\x n ) 
specified by the M-H algorithm, namely if 



a(x'\x n ) = rain j-i[x 2 (a;') - X 2 (x n )}, 1 j > u 



(B2) 



where u is a random number generated from a uniform distribution between and 1. For the Gibbs sampler, 
the trial state x 1 is generated by substituting only a subset of parameters from the state x. In particular, 
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since Rt, S t and 7 are correlated (Eq. [9]), we choose to change these three parameters at the same time, 
while varying the inclination and the position angle independently. At each step, the parameter(s) to be 
updated are then randomly modified with a transition probability of 

= exp 

where the variance (3^ defines the variability interval for each parameter. The model fitting results depend 
neither on the choice of the transition probability q(x'\x) nor on how the Gibbs sampler is implemented. 
However these strongly effect the efficiency of the model fitting and must be carefully chosen to allow a 
fast convergence toward the equilibrium distribution. Here, each parameter is allowed to vary over a large 
interval and the variance [3^ is chosen so that the overall acceptance rate is close to the optimal value of 
-0.25 (Gelman et al. 2003). 

To select the initial state of the chain we adopt two different approaches. First, the initial parameters for 
the disk are set based on the Gaussian fit to the observations (Table |3|) ; the initial inclination and PA are 
calculated from the aspect ratio of the emission, the transition radius Rt is assumed to be the half width 
along the major axis at half maximum, p is set to be 0.5 and £ t is randomly chosen between 0.1 and 1000 
g/cm 2 . This choice is usually close to what provided by the best fit model and enables a quick convergence 
of the chain. Unfortunately like most algorithms developed to minimize the x 2 > the MCMC can be trapped 
in local minima if they are separated by sufficiently high \ 2 barriers. To compensate, additional randomly 
initialized MCMC chains were run to verify that they all converge to the same final state. 

Once the equilibrium distribution of the MCMC chain is sufficiently well sampled - usually requiring a run 
of a minimum of 10 4 models - the distribution of each parameter is obtained through marginalization, i.e. by 
integrating the posteriori distribution (now equivalent to the MCMC equilibrium distribution) over all the 
parameters except the one in which we are interested. We show in Figure [18] the obtained posteriori distribu- 
tions for the model parameters. These are generally consistent with a normal profile, supporting our decision 
to adopt a Gaussian transition probability q(x'\x), and enables expressing the parameter uncertainties in 
terms of the standard deviation a of the probability distribution. 



2/32 



(B3) 



Cumulative and total disk mass 



We integrate the surface density (Eq. [9]) to derive the cumulative disk mass, i.e., the disk mass contained 
within a radius R at time t. Since Ri n <C Rt we obtain 



M d (R,t) = 47rS t i? 2 e 1/2(2 " 7) x \ 1 - exp 



1 / R 

"2(2-7) 



(2- 7 ) 



The total disk mass at the time t is therefore obtained for R — > 00 in the form 

M d (t) = 4 7 rS tJ R t 2 e 1 / 2 ( 2 -T). 



(CI) 



(C2) 



Using the last two equations we derive that the disk mass contained inside the transition radius Rt is 39% 
of the total disk mass if 7 = 1 and 22% if 7 = 0. 

Using Eq. IC1I and Eq. IC2I we can also demonstrate that R\ is the radius containing ~63% of the initial 
disk mass Md(t = 0). We can write 



M d (fli,0) 
M d (0) 



1 — exp 



1 



2(2-7) \Rt,o 



Ri 



(2-7) 



(C3) 
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where the transitional radius at t = 0, given by Eq. [8l is 



Rt.o = Ri 



1 1/(2-7) 



.2(2-7) 

Substituting this latter equations in Eq. IC3I we obtain 

M d (R u 0) 



M rf (0) 



1 -e- 1 ~ 0.63 



We finally calculate that about 90% of the initial disk mass is contained within 2Ri . From Eq. 
Eq. IC4I we can write 

M d (2R 1 ,0) _ i 2 (»-t) 
M d (0) 6 
which is equal to 0.86 for 7 = 1 and 0.98 for 7 = 0. 



(C4) 

(C5) 
and 
(C6) 



D. Mass accretion rate on the central star 

From Eq. O and [H] the surface density S t at the transition radius Rt can be expressed as: 

1 



Et = JL T -{V3-7)/(2-tO f #-V x exp 
3iTisi \R t I 



2(2-7) 



Writing the initial disk mass Md(0) in the form 



MM = l° Rl 



3 v x 2 - 7 

we can eliminate the ratio C/z/i from Eq. IDll to obtain 



= -Mrf(0) 1/2(2 _ 7 ) 1/2(2-7) 

^Rl 



Substituting the definition of R t from Eq. [5J S t takes the form 



S 4 oc 



Mj(0) y_ 5 /(2(2--y)) 



Starting from the expression of the mass flow (jHartmann et al. Ill998j ) 



. (5/2-7) 

M(r, t) = CT~ (=-7) exp 



,(2-7) 



2(2- 7 )r( 2 "^) 



(Dl) 



(D2) 



(D3) 



(D4) 



(D5) 



and substituting the expression for the constant C derived from Eq. ID2I we can express the mass flow as 



M(R,t) = 



M d (0) 



1-1 * 

Rt 



(2- 7 )' 



2(2-7)* 

where t s is the viscous time scale defined by the Eq. [6] and 



jT- (5/2- 7 )/(2-7) 



exp 



i?(2-7) 



2(2 - 7 )i? t 2 - 7 



R R Rt R 
Ri Rt Ri Rt 



T 



1/(2-7) 



.2(2-7). 



(D6) 



(D7) 
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Note that for 7 = 1, Eg. ID6I reduces to Eq. 35 in Hartmann et al. ( 19981) . 

At radii smaller than the transitional radius Rt, the material within the disk moves inward and finally 
falls onto the central star. The mass accretion rate M acc (t) is given by Eq. ID6I with R equal to the radius 
at which the disk is truncated by the accretion process. Since this radius is of the order of a fraction of AU 
and much smaller than the transitional radius R t , M acc takes the form 



M acc (t) 



MM 
2(2 - 7 )t 



_ T _(5/2- 7 )/(2- 7 )_ 



(D8) 



E. Derivation of a in the case 7 7^ 1 

In this section we derive the value of a for the similarity solution of the disk surface density in the general 
case 7^1. We adopt the classical a parameterization of the disk viscosity in the form (Shakura & Sunyaev 
1973) 

v = ac s H. (El) 

If the disk is in keplerian rotation, vertically isothermal and in hydrostatic equilibrium (see the discussion 
in § 15. 2j) , the sound speed is given by 



c s = H ■ Vt = H ■ y/GM+fR 3 (E2) 

with 

where kf, is the Boltzman's constant, [i — 2.33 g mol -1 (Ruden and Pollack 1991) is the mean molecular 
weight of the circumstellar material, mu is the proton mass, G is the gravitational constant and Ti(R) is 
the disk interior temperature. Substituting Eq. IE2l and lE3l in Eq. IE1[ we can express the disk viscosity as a 
function of the radius in the form 

v{R) = k ±== ■ aR 3 / 2 Ti(R). (E4) 

M m £T V GM* 

The comparison of Eq. \W\ and Eq. IE4I provides therefore an expression of a in the form 

3xio- 2 / t. V7 Ri V 2 ~ 7) ( R \ (7 - 3/2 V 



< R ^^W\jMVr) UxuJ UauJ UaSp Ti {R) - (E5) 



'0 
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Table 1 
Sample properties 



Object 


a(2000) 


(5(2000) 


ST 


L* 


T* 


log(L acc ) 
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R* 


log(Macc) 
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3 00 

o.uu 


/LOfiO 
^UOU 


A A 


Q P. *7 
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ft fi8 
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U. 1 
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04:43:03.09 


25:20:18.59 


M0 


0.28 
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-0.98 


1 
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-8.33 


0.6 


3.0 
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04:39:17.78 


22:21:03.52 


K5 


0.74 


4350 


-1.75 


1 


1.60 


-9.17 


0.7 


1.8 


RY Tau 


04:21:57.41 


28:26:35.56 


Kl 


7.60 


5080 


0.20 


5 


2.92 


-7.11 


2.0 


0.5 


UZ Tau E 


04:32:43.07 


25:52:31.14 


Ml 


0.90 


3720 


0.17 


5 


2.28 


-6.90 


0.3 


0.4 


GM Aur 


04:55:10.98 


30:21:59.38 


K7 


0.74 


4060 


-1.15 


1 


1.4 


-8.55 


0.5 


1.0 


GSS 39 


16:26:45.00 


-24:23:07.70 


Ml 


1.20 


3720 


-0.40 


2 


2.64 


-7.43 


0.3 


0.1 


SR 24 S 


16:26:58.50 


-24:45:36.90 
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2.50 


4170 


0.03 
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3.03 


-7.13 


0.4 
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TW Hya 


11:01:51.91 


-34:42:17.02 
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-21:57:21.88 


Al 


36.0 


9500 


0.40 


4 


2.2 


-7.12 


2.3 


5.0 



Note. — In column (5) we report the stellar luminosity in solar luminosities, in column (6) the stellar temperature 
in K, in column (7) the accretion luminosity in Lq/yt, in column (9) the stellar radius in solar radii, in column (10) 
the mass accretion rate in M Q /yr, in column (11) the stellar mass in solar masses and in column (12) the stellar age 
in Myr. 



References. — I'll lHartmann et al. I (|l998l ). (2) iNatta et all (|2006l ). (3) iMuzerolle et al. I (|2000l ). (4) 
Garcia Lopez et al. I ((2006J), (5) Kenyon and Hartmann (1995), (6) Calvet & Gullbring (1998), (7) Muzerolle et 
al. (1998) 



Table 2 

Summary of CARMA continuum observations at 230 GHz 



Object 


Config. 


Beam FWHM (") 


Beam PA (°) 


Date 


Phase Calibrators 


CY Tau 


C 


1.05x0.72 


77 


2007 Nov 12 


3C111, 3C123 


DG Tau 


C 


0.87x0.78 


-50 


2007 Sep 30 


3C111, 0530+135, 3C123 




c 


0.83x0.64 


-72 


2007 Oct 08 


3C111, 0530+135, 3C123 




B 


0.43x0.27 


-74 


2007 Dec 14 


3Clll, 0530+135, 0510+180 


DM Tau 


C 


0.82x0.60 


-78 


2007 Nov 05 


3C111, 0530+135, 3C120 


DN Tau 


C 


0.80x0.58 


-76 


2007 Nov 05 


3C111, 0530+135, 3C120 


DR Tau 


C 


0.92x0.76 


-83 


2007 Oct 09 


3C111, 0530+135, 0449+113 




C 


0.84x0.70 


82 


2007 Oct 24 


0530+135, 0449+113 




B 


0.46x0.34 


39 


2007 Feb 06 


0530+135, 0449+113 


GO Tau 


C 


0.87x0.65 


88 


2007 Nov 07 


0530+135, 3C123 


LkCal5 


c 


0.83x0.70 


70 


2007 Oct 27 


3C111, 0530+135, 3C123 


RY Tau 


c 


0.89x0.74 


-61 


2007 Oct 01 


3C111, 3C123 




c 


1.14x0.60 


-73 


2007 Oct 22 


3C111, 0530+135, 3C123 


UZ Tau E a 


c 


0.82x0.69 


79 


2007 Oct 27 


3C111, 0530+135, 3C123 


GSS 39 


c 


1.42x0.85 


-6 


2008 Apr 12 


1625-254, 1733-130 


SR 24 S b 


c 


1.45x0.91 


-6 


2008 Apr 13 


1625-254, 1733-130 



a UZTau W, the other component of the UZTau system, was detected at 4cr level with an integrated 1.3 mm 
flux of about 30 mJy (see Fig. QJ 

The other component of the SR 24 binary system, SR 24 N ([Patience et al. 2008), was not detected. 



Table 3 

Properties of the 1.3 mm dust emission 



Object 


Flux (mJy) 


Source size 




PA (°) 


R g (AU) 






FWHM (") 








CY Tau 


117+20 


1.25x0.60 


61 


150 


230 


DG Tau 


317+28 


0.52x0.46 


28 


11 


95 


DM Tau 


90+8 


0.89x0.82 


22 


20 


160 


DN Tau 


93+8 


0.68x0.53 


39 


86 


125 


DR Tau 


109+11 


0.48x0.39 


36 


108 


90 


GO Tau 


57+8 


0.86x0.69 


37 


107 


160 


LkCal5 


119+15 


1.39x0.71 


59 


55 


250 


RY Tau 


227+20 


0.63x0.38 


53 


23 


115 


UZ Tau E 


126+12 


0.88x0.64 


43 


66 


160 


GM Aur 


189+15 


1.38x1.08 


45 


106 


270 


GSS 39 


282+20 


1.33x0.87 


49 


116 


260 


SR 24 S 


197+17 


1.64x0.96 


54 


60 


320 


TW Hya 


543+45 


1.06x1.04 


12 


89 


76 


MWC 275 


705+12 


1.49x1.12 


41 


135 


250 



24 



Table 4 
Dust properties 









Umax 


=100mm 


(7= 


=3 


Object 


a 


P 


q 


fci.3 a 


a b 


fcl.3 a 


CY Tau 


2.6 


0.7 


3.5 


0.6 


1 


3.9 


DG Tau 


2.3 


0.5 


3.1 


0.2 


22 


0.39 


DM Tau 


2.9 


1.1 


4.0 


1.8 


0.16 


6.6 


DN Tau 


2.7 


0.8 


3.7 


1.2 


0.55 


5.6 


DR Tau 


2.4 


0.6 


3.3 


0.3 


1.9 


2.6 


GO Tau 


3.4 


1.5 


4.5 


0.4 


0.13 


6.8 


LkCal5 


3.5 


1.7 


4.7 


0.3 


0.13 


6.8 


RY Tau 


2.5 


0.7 


3.5 


0.6 


1.0 


3.9 


UZ Tau E 


2.6 


0.7 


3.5 


0.6 


1.0 


3.9 


GM Aur 


3.1 


1.3 


4.2 


1.1 


0.14 


6.7 


GSS 39 


2.8 


0.9 


3.8 


1.6 


0.3 


7.4 


SR 24 


2.6 


1.1 


4.0 


1.8 


0.16 


6.6 


TW Hya 


2.5 


0.8 


3.6 


0.9 


0.55 


5.6 


MWC 275 


2.9 


1.0 


3.9 


1.9 


0.185 


7.0 



a 1.3 mm dust opacity is in cm 2 per gram of dust. 
b dust size is in mm. 



Table 5 
Model fitting results 



Object 
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Rt 


s t 
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Log(M d ) 


Rd 
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(deg) 


(AU) 


(g/cm 2 ) 




(M e ) 
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Fig. 1. — : 1.3 mm dust continuum images of thc^fjisks observed with CARMA. Contours start at the 
significance levels given in each panel and are separated by that same amount. The exception is UZ Tau E/W 
where a cross indicates the position of UZ Tau W and the dotted contour corresponds to the 4tr level. Beam 
sizes and PA are listed in Table O Integrated fluxes and source sizes are given in Table [3] 
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Fig. 2.—: 1.3 m m dust continuum maps of GM Aur, TW Hya and MWC 275 (from iHughes et al. 



Isella et al. 1 120071) . For GM Aur and MWC 275 contours begin at the 3cr level and are separated by the 
same amount. For TW Hya contours start at and are separated by the 10a level. Beam sizes and PA are 
listed in Table [2] Integrated fluxes and source sizes are given in Table [3] 
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Fig. 4. — : Disk surface density as denned by Eq. [5]for 7 = -1, and 1. The transition radius and the 

normalization are fixed at Rt — 30 and Et = 50 g/cm 2 . 
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Fig. 6. — : Flux density distribution of the observed sample between 0.4 and 8 mm. Filled squares corre- 
spond to the CARMA observations discussed in this paper, and open squares correspond to measurements 



spond to the LAKMA observations discusse d m this paper, ana open squares c orrespond, to mea surement 
available in literature jDutrev et al. 1996; iKitamura et al. 1 120021 I Andrews et al. 1 12005 1: iRodmann et al 



2006; llsella et al. 1 120071 ). The solid curves correspond to disk models that best fit the spatially resolved 
observations at 1.3 mm (see Table [SJ. For each object the values of a and of the dust opacity slope (3 are 
indicated in the relative panel. 
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Fig. 7. — : Radial profiles of the disk surface density, cumulative flux at 1.3 mm, cumulative mass and disk 
temperature for the best fit model from Table The shaded region in the surface density panels corresponds 
to the la uncertainties obtained from MCMC fitting. The vertical dashed lines indicate the spatial resolution 
of the observ ations. In the p anels of column 3, the horizontal dashed line at 0.02 Mq is the minimum mass 
solar nebula (jHavashi 1 1 1 9 8 lh . and the disk cumulative mass is expressed in units of 10 -2 Mq. The disk 
interior temperature Tj is represented by a thick line in the panels of column 4 while the disk surface layer 
temperature is a thin line. 
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Fig. 7. — : continued 
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Fig. 7. — : continued 
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Fig. 8. — : Comparison between the observed correlated flux (dots) and the best fit model prediction (solid 
line) as a function of the bcsclinc lcnght deprojected using disk inclinations and position angles listed in 
Tab. [5] The histogram in the lower part of each panel shows the expected signal in case of zero flux. 
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Fig. 9. — : Maps of the residuals calculated by subtracting the best fit model from the observations. The 
contours start at the 2a level and are separated by la. Cross indicates the position of the source and the 
orientation of the disk. The smaller cross in the UZ Tau E panel indicates the position of UZ Tau W. 
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Fig. 9. — : continued 
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Fig. 10. — : The disk transition radius R t as a function of the stellar age computed using DM97 models. 
Filled squares show the 10 disks located in Taurus- Auriga star forming region. For this sub sample, the 
correlation coefficient between R t and the stellar age is 0.98 and the probability that the data are randomly 
distributed is less than 0.1%. The solid line corresponds to R t = R +C-t' n with r) = 0.5±0.4, R a = 17±10 AU 
and C = 37 ± 20. 




Fig. 11. — : The mass accretion rate as a function of stellar age computed from DM97 models. Filled squares 
show the 10 disks located in Taurus-Auriga star forming region. The correlation coefficient between M acc 
and the stellar age is -0.62 and the probability that the data are randomly distributed is about 2%. The 
solid line corresponds to M acc oc t~ 1A . 
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Fig. 13. — : Radial profile of the gravitational instability parameters Q. The solid line corresponds to 
DG Tau, the long dashed line to LkCa 15 and the short dashed lines to other sources in the sample. The 
grey region (Q < 1.5) indicates where the disk is gravitationally unstable. 
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Fig. 14. — : Derived values of the stress parameter a, for the observed circumstellar disks grouped according 
to 7 > (upper panel), 7^0 (middle panel) and 7 < (lower panel). 
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Fig. 1 5. — : Spectral ener g y dist ri bution of LkCal5 . Flux measurements from 2MA SS, lKenvon fc Hartmann 



(|l995l ). iKitamura et al. I (|2002t ). iRodmann et al. I (|200d ). lAndrews et"al~l (|2007t) are represented by open 
squares, open circles and open triangles respectively. Our CARMA data are shown as filled circle, while data 
between 5 and 14 [im are from the Spitzer IRS archive. The solid line is the SED for the disk model that 
fits our CARMA 1.3 mm continuum observations (Tab. [5]). This comprises the stellar photosphere (dotted 
line), a "puffed-up" inner rim (long-short dashed line), a disk surface layer (short dashed line) and a disk 
midplane (short dashed line). 
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Fig. 16. — : Position on the H-R diagram of the observed sources. The upper panel show the theoretical 
models by D'Antona & Mazzitelli (1997) while the lower panel the models of Baraffe et al. (1998). The 
dashed lines correspond to the stellar isocrones for ages in Myr as labeled at the left end of the lines while 
the solid lines correspond to the stellar evolution sequence for stellar masses between 0.2 and 1.7 M Q as 
labeled at the lower end of the lines. 
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Fig. 17. — : Upper panel: stellar ages derived from DM97 versus the ages derived from the B98 models. The 
solid lines correspond to equal ages. The comparison between stellar masses is shows in the lower panel. The 
open squares identify the dynamical stellar masses derived by Simon et al. (2000) for DM Tau, CY Tau, 
GM Aur and LkCa 15 
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Fig. 18. — : Normalized probability distributions for the disk inclination (°), PA (°), transition radius Rt 
(AU), surface density S t (g/cm 2 ) and 7 obtained from the MCMC fitting process as discussed in Appendix IB1 
The solid line corresponds to the best fit Gaussian distribution used to derive the parameter uncertainties 
reported in Tab. [5] 
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